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Abstract 

In recent years modelling crowd and evacuation dynamics has become very important, with increasing 
huge numbers of people gathering around the world for many reasons and events. The fact that our global 
population grows dramatically every year and the current public transport systems are able to transport 
large amounts of people, heightens the risk of crowd panic or crush. Pedestrian models are based on 
macroscopic or microscopic behaviour. In this paper, we are interested in developing models that can be 
used for evacuation control strategies. This model will be based on microscopic pedestrian simulation 
models, and its evolution and design requires a lot of information and data. The people stream will be 
simulated, based on mathematical models derived from empirical data about pedestrian flows. This model 
is developed from image data bases, so called empirical data, taken from a video camera or data obtained 
using human detectors. We consider the individuals as autonomous particles interacting through social 
and physical forces, which is one approach that has been used to simulate crowd behaviour. 

Typical questions in pedestrian motion planning and design are: How many minutes would be needed 
to evacuate a football stadium or a pop concert arena? Where is the best place for pedestrian information or 
a poster wall to be placed? How many exit signs should be distributed over the area and where shoidd they 
be placed? How many minutes does a pedestrian need to go from point A to point B in the overcrowded 
places? Are transfer times still realistic when the pedestrians are not familiar with the layout?. For better 
planning in this area a greater understanding of human behaviour is required and one way to achieve 
this is to improve the tools available for pedestrian planning (such as pedestrian micro simulation). In 
particular, particle-based simulation has become an attractive paradigm for modelling and simulating 
pedestrian decision-making and behaviour because it supports a one-to-one correspondence between the 
subject of observation and the simulated particle. 

The target of this work is to describe a comprehensive approach to model a huge number of pedestrians 
and to simulate high density crowd behaviour in overcrowding places, e.g. sport, concert and pilgrimage 
places, and to assist engineering in the resolution of complicated problems through integrating a number 
of models from different research domains. 

Keywords: Pedestrian Dynamics - Crowd Simulation and Modelling, Crowd management and 
pedestrian safety. Crowd control. Objects tracking. High Density Pedestrian Flow. 


1 Introduction 


I n recent years, several models for the move¬ 
ment of crowds have been proposed. One 
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can distinguish between two general ap¬ 
proaches: microscopic and macroscopic mod¬ 
els. In the microscopic scope (range), pedestri¬ 
ans are treated as individual entities (particles). 
The particle interactions in this model are de¬ 
termined by physical and social rules as well 
as the interactions between the particles with 
their physical surrounding. In this context we 
cite, the social-force models (see m and the 
references therein), the cellular automata, e.g 
mist queuing models e.g. fU, or continuum 
dynamic approaches like f5]|. Tor an extensive 
survey of different microscopic approaches we 
refer to 0. 

The origin of the microscopic models for 
pedestrian behaviour can be traced back to the 
work of Reynolds in the years 1987 and 1999 
00. This model was extended by the socio¬ 
logical factors introduced in (9j, psychological 
effects in ffTOl , social forces IHT1IT2 , 13,141, and 
other models of pedestrian behaviour [15, 16, 
[17]. Many methods have been developed to 
describe the short range forces responsible for 
the collision avoidance between nearby pedes¬ 
trians. These include geometrically-based al¬ 
gorithms fll8l Il9l l20l |2T1 l22l l23ll . grid-based 
methods li24l [25] [26l| and force-based methods 


In contrast to microscopic models, macro¬ 
scopic models treat the pedestrians flow as con¬ 
tinuum homogeneous mass that behaves like a 
fluid without considering the behaviour of sin¬ 
gle individuals. Many methods use concepts 
familiar from fluid and gas dynamics, see Ii30l . 
Modelling pedestrian dynamics as a fluid can 
be traced back to 1131] |32l . More recent models 
are based on optimal transportation methods 
ll33l . mean field games (see If34l for a general 
introduction) or non-linear conservation laws 


2 Pedestrian forces and mo¬ 
tion (in PedFlow) 

PedFlow is a microscopic simulation model 
developed by Lohner Simulation Technologies 
International, Inc. (LSTI) 1361 . In this model 


individuals are treated as self-driven particles 
that interact through social and physical rules. 
Each person in the pedestrian flow has a de¬ 
sired position and desired velocity and adapts 
his or her current velocity according to the sur¬ 
rounding neighbour; furthermore, each crowd 
member simultaneously tries to avoid collision 
with other crowd members and any environ¬ 
mental boundaries. The current velocity and 
direction of a single individual in the crowd is 
a result of the circumstances and social inter¬ 
actions. These social forces can be influenced 
by the environment, other people and internal 
states. 

Force plays an essential role in every crowd 
model. It can be classified in two major cate¬ 
gories, internal and external forces. Forces are 
characteristics of crowds and crowd disasters. 
Force has a direct effect on movement: for ex¬ 
ample force drives people out of rest to their 
desired velocities, or force can make deviate 
people from their desired direction. In dan¬ 
gerous crowd situations, force can be the key 
to decode the fatal consequences (e.g. injuries 
and death). Many force effects are summarized 
and documented by Fruin (37|. 

In this approach the pedestrians are treated 
as "particles'' moving according to Newton's 
law: 

dv 7 dr _ 

-.=° in 


dt 


Here nip, v and r denote the mass, velocity and 
position of the pedestrian, and / the sum of all 
forces exerted by it or acting on it. The basic 
unknown that requires intensive modelling ef¬ 
forts in these equation is the force-vector. In 
the sequel, we will discuss the main forces that 
can be identified acting on a pedestrian. 

2.1 PedFlow force models 

In the PedFlow model, the forces that acceler¬ 
ate a pedestrian are the result of internal (will) 
and external (collision, signals, etc.) forces. An 
"individual alert" tries to avoid collisions before 
they happen, resulting in vanishing external 
forces. However, if the pedestrian density in¬ 
creases, collision forces will appear, and the 
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ability to move freely will be impaired. The 
difference between external and internal forces 
is too rigorous and only serves descriptive pur¬ 
poses in the present context. So, what are the 
forces that act on the pedestrian? 

• Internal forces (active, intentional) 

- Will force (get there in time) 

- Pedestrian collision avoidance 
forces 

- Wall/Obstacle collision avoidance 
forces 

• External forces (passive, suffered) 

- Physical pedestrians collision forces 
(physical contact) 

- Physical wall/obstacle collision 
forces (physical contact) 

2.2 General description of will force 

Each particle is programmed to have a final 
goal that it wants to reach. The motivation to 
move towards this goal is given by the driving 
force, or so called will force f w m and denotes 
the self-driven force that drives an individual 
towards their desired velocity (see fig. [lj. The 
self-driven force is based on a simple error 
correction term consisting of the difference be¬ 
tween an individuals desired velocity, and a 
desired direction z , and their actual velocity. 
This difference in velocity is corrected over a 
specified time interval T, referred to as the re¬ 
laxation time, that corresponds to the finite 
amount of time that is required for people to 
react and physically change their velocity. 

In the following, we denote by will force 
the force that will accelerate (or decelerate) a 
pedestrian to achieve the velocity it desires. 
Given a desired velocity v d and the current 
velocity v, this force will be of the form 

fwill=gw{vd-v), (2) 

where g w — f(r) is a function of "list of cir¬ 
cumstances" 



Figure 1: The will force f w m denotes the self-driven force 
that drives an individual towards their desired 
velocity. 


The modeling aspect is included in the func¬ 
tion g w , which in the non-linear case may itself 
be a function of v d — v. Suppose g w is constant, 
and that only the will force is acting. Further¬ 
more, consider a pedestrian at rest. The pedes¬ 
trian could start from rest Vq — 0 and wants 
to reach the desired velocity v d . In this case, 
we have: f w m = gw{v d — v ), with u(0) = 0 we 
have m p f t = f wlU = g w (v d - v). The solution 
of this equation is given by v(t) — v d (l — e~ at ) 
with ^(0) = v d tx — v d /T and a = g w /m p — 
1 / r We can thus obtain g w by measuring the 
time required to reach a percentage (e.g. 90 
percent) of the desired velocity, starting from 
rest (see fig. |2j. This time is typically in the 
range of 0.5 - 1.0 s, but obviously depends 
on the current state of fitness, stress, climate 
and terrain condition, and desire to reach a 
goal. We can define the function g w — m p /t 
via "relaxation time" 0.5 < r < 1.0 s. Typical 
values for the desired velocity v d and relax¬ 
ation time r are v d = 1.35 m/s and r = 0(0.5 to 
1.0 s). Medicinal facts and human experience 
show that older pedestrians walk more slowly 
than younger pedestrians. This effect has a 
significant impact on the "relaxation time" of 
a single pedestrian, hence it is easier to assess 
that younger pedestrians reach their desired 
velocity faster than older pedestrians. 
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Figure 2: Speed-up of an isolated particle to the desired 
velocity. This difference in velocity is corrected 
over a specified time interval r , which is re¬ 
ferred to as the relaxation time, that corre¬ 
sponds to the finite amount of time that is re¬ 
quired for people to react and physically change 
their velocity. 

2.3 Internal forces 

2.3.1 Internal forces -far range 

In normal cases, the will force f w m takes the 
direction of the desired velocity as shown in 
figure [3] 



Figure 3: Desired velocity. In the real world, each pedes¬ 
trian has a desired velocity iq towards a de¬ 
sired destination and adapts the current veloc¬ 
ity according to the environmental situation. 

Traditionally the will forces are specified 
by a desired velocity of the pedestrian and the 
trajectory. The will force will be implemented 
as a vector on the simulation model. However, 
in the motivation to reach certain places, there 
are tasks in which the desired motions of indi¬ 
viduals are specified by the state of system (list 
of circumstances) rather than time. For such 
tasks, desired velocity field control has been 


proposed. In this chapter, a method to extend 
our information about the forces that drive the 
individuals from rest to their desired velocities 
is evaluated. 

• In the case of group movement, the will 
force takes the direction of the group 
leader. 


f will = will force 



Figure 4: "Follow me". 

Figure |4] illustrates a group of pedestri¬ 
ans following their leader. A "group” is 
defined as a physical collection of people 
following the same route, but who may 
or may not be part of the same social 
group, and a "subgroup" is defined as 
people within the same physical "group" 
who want to stay together, like friends 
or family members. Several studies have 
revealed that smaller subgroups consti¬ 
tute the majority of the people in a crowd. 
But very few studies are able to model 
the "subgroup" behaviour. A particu¬ 
lar subgroup concerns pedestrians who 
hold hands. It is recognized that holding 
hands is the most effective way of keep¬ 
ing children safe from traffic injury [38]. 
All the group members take the desired 
position and the desired velocity of the 
group leader. 

• In the case of an emergency or dangerous 
situation, the will force takes the direc¬ 
tion of the nearest exit and this has dra¬ 
matic consequences on the behaviour of 
pedestrians at the exit. They block each 
other and this builds up a rigid arc from 
which no one can escape (see fig. [5] and 

0 - 
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Figure 5: Go to nearest exit. 


Figure [5] demonstrate the situation of pedes¬ 
trians in case of emergency or dangerous sit¬ 
uation. The will force points to the direction 
of the nearest exit. This picture illustrates the 
situation when people take the proper precau¬ 
tion to prevent potential hazards when they are 
alerted. For example, in case of an emergency, 
evacuating the building immediately by using 
the closest exit and going to the designated 
meeting location. Panic breaks out, if many in¬ 
dividuals attempt to escape an emergency area 
at the same time, often the emergency exits 
are blocked in a short time, and the situations 
have a fatal consequence for survival and leads 
to injuries or death through crushing people 
in a crowd. Real world example (fire safety: 
General Evacuation Procedures): Learn the lo¬ 
cation of the two nearest exits from your work 
area (often these are down staircases) ll39| . 



Exit 


Figure 6: Scheme of an arc formation at openings and 
bottleneck areas. The full ellipses represent 
pedestrians mooing in the same directions. The 
pedestrians block one another in the direction 
of the movement and congest the passage area. 

The priority in the emergency or dangerous 
situations like fire in a building or earthquake 
are designed to reach each particle the nearest 
exit in a short time. In first step the particle 
choose their priority destination after that they 
can find the beeline towards the nearest exit. 

2.3.2 Internal forces -intermediate range 

(a) Collision avoidance forces 
In the collision avoidance process, autonomous 
particles need to discover the environment to 
avoid static and dynamic obstacles. At any 
time the distance between each obstacle, wall, 
and the particle i are computed, go both of 
them closer to each other, then calculate the an¬ 
gle between particle i's desired direction and 
the line joining the center of particle i and the 
obstacle. This effect can be seen within the 
rectangle of influence (figure [7|. The distance 
and the angle play a significant role to assess 
how relevant the obstacle is to the trajectory. 
When travelling the environment, particles also 
update their perceived crowd density, which is 
necessary for their decision making process. 

We denote by intermediate range forces that 
change the motion of an individual in order to 
avoid an encounter with another. The obser¬ 
vation employed here is that all of us will try 
to avoid a collision by moving away from an 
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encounter long before it happens. The nature 
of these forces is such that they act mainly in 
a direction normal to the current motion of 
the individual, and therefore do not lead to a 
significant decrease in velocity. Moreover, this 
force is only active if a collision is sensed. 



Figure 7: Collision avoidance rectangle of influence. Il¬ 
lustration of pedestrian-pedestrian collisions 
avoidance forces. This force is referred as re¬ 
pulsive force and reaches its max. value at the 
shortest distance between the particles. 

Using as a starting point the situation de¬ 
picted in figure |7j given two pedestrians with 
current coordinates and velocities f\, V\ and ? 2 , 
z> 2 , an encounter may be computed by eval¬ 
uating the time increment At at which the 
distance between the two is minimized, i.e., 
[?i + A tv\ — (?2 + Afz^)] 2 —> min. This results 
in [(ri - r 2 ) + At m (fi ~ v 2 )\{vi - Vi) = 0, or 
A t m — — . The minimum distance 

is given by S min = |(r a - f 2 ) + At m (v x - v 2 )\ . 
Obviously, the force will only become active 
if At > 0. The general nature of this force is 
such that it will decrease with distance. How 
exactly this decrease function looks like is un¬ 
known at the present time (it may even be 
random). We have used a function based on 
the normalized distance between the pedes¬ 
trians d — rz ^ , where R\ is the charac¬ 
teristic radius of pedestrian 1, which is of the 
form: / = f max ^ , ■ This simple repulsion 

force may be refined further by considering 
the directions along and normal to the current 
velocity vector v\, denoted by e t and e n , re¬ 


spectively. We can then split the normalized 
distance d into a tangential and normal com¬ 
ponent: dt — ihhvzlhl' d n — !_ This 

leads to tangential and normal forces of the 
form: /= -f max The 

value of f max can be related to the relaxation 
time, i.e. it is not too difficult to obtain. We 
have used f max /mp = 0(4) « 4 [m/sec 2 ]. 

(b) Avoiding a wall 

Avoidance forces /repulsive are calculated only 
for relevant obstacles, walls and agents: i.e. 
those falling within the rectangle of influence. 
The avoidance force for obstacle W is illus¬ 
trate in figure |8j The calculation is done in 
the same way as for the force of the pedestrian- 
pedestrian repulsion and the same formulas 
are used. Only the pedestrian-pedestrian dis¬ 
tance has to be replaced by the distance of the 
pedestrian from the wall. 



pedestrian 

trajectory 


Figure 8: Avoiding a wall. This picture illustrate the 
pedestrian-obstacle collision avoidance. The 
blue solid arrow shows the direction of the re¬ 
pulsive acting force f repulsive a) id the dashed 
arrow indicates the trajectory path of the parti¬ 
cle. 

2.4 External forces 

2.4.1 External forces -near range 

• Pedestrian collision: For more realistic 
crowd simulation and for a description 
of realistic counter-flows and overtaking 
behaviours, many attributes influencing 
the behaviour of people in the real world 
will be considered. This approach has the 
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ability to simulate human behaviour by 
setting parameters obtained from many 
observations related to real human move¬ 
ment. In avoiding an obstacle the factors 
affecting the tangential forces are: 


- The distance between the particle 
and any obstacles. 


- The relative direction and move¬ 
ment of the agents i.e. the direction 
between the desired velocities vec¬ 
tors (z>| and z >2 for example). 


- Crowd density: If an agent enter 
the interaction rectangle, a tangen¬ 
tial force (see fig. [7) will be ac¬ 
tive to change the direction of move¬ 
ment and make a curve in the path 
to avoid collision. The angle be¬ 
tween two agents is used to sim¬ 
ulate the individual decision mak¬ 
ing how to react to a previous col¬ 
lision. For example, in a walking 
floor, if there is enough space be¬ 
tween us and another person, none 
of us would change the direction. 
In case of many persons walking in 
the corridor, the majority of people 
tend to move to the right. Therefore, 
when the velocity vectors are almost 
collinear, the tangential forces point¬ 
ing to the right. In case an agent 1 
detects agent 2 and agent as possible 
obstacles, then the system computes 
the distance vector towards agent 1 
as follow (d 2 \). Agent 2 is farther 
away from 1, but since it is moving 
against agent 1, the system gives this 
obstacle higher priority. 



Figure 9: Pedestrian collision. Schema illustrate 
pedestrian-pedestrian body contact or colli¬ 
sions forces. This situation indicate a body 
friction or body deformation. 


Once a pedestrian has moved close 
enough to obstacles or other pedestrians, 
his velocity will be decreased markedly, 
(see fig. |9j. Unlike the long-range forces, 
these forces act in the direction of the 
normalised difference vector z — r f 1 ■ 
We have used, as before, a force of the 
form f — —fmax (i+ d 2 ) Z, although the ex¬ 
act nature of this force is unknown at the 
present time. As before, the value of f max 
can be related to the relaxation time. We 
have used f max / m P = 0(4) ~ 4 [m/sec 2 ]. 


• Collision with a wall or obstacle: A 
stream of pedestrians will be enclosed 
by the walls in or around the structures 
in which it occurs (see fig. 10 1. We have 
assumed that this force acts in the di¬ 
rection of the gradient of the distance- 
to-wall function d w (x). This function, 
which measures the closest distance to 
a wall from any given location r, is as¬ 
sumed to be known. In practice, it is 
evaluated in a pre-processing step. As be¬ 
fore, we have used, for lack of any deeper 
knowledge, a function of the form: f w — 
Jmax (i_j -(rf-p/R) 2 ) ' dzvr with fmax/ tit p — 

0(4) ~ 4 [m/sec 2 ]. 
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contact area 



Figure 10: Repulsive contact force. The solid blue arrow 
indicates the direction of the contact repulsive 
force. 

f denotes the individual-individual inter¬ 
action force or individual-obstacle force, 
f w denotes the collision force due to the 
wall (out of a total of M wall surfaces). 

2.5 Data structures 

2.5.1 Data carried by the individual 

The data concerning a single individual are 
classified in three categories: physical, emo¬ 
tional and ethnic data. This data, which in¬ 
cludes items such as the pedestrian's height, 
width, fitness, nationality, current physical and 
mental state, as well as the current motion 
data (location, velocity, etc.) is stored in ar¬ 
rays attached directly to the pedestrians. In 
this way, the data becomes easily accessible 
at any time. In order to save storage, the in¬ 
dividual's height, width and fitness are taken 
from a data-base, i.e. through a table look-up. 
In this table, the data belonging to a limited 
number of representative groups is stored, and 
the individual is tagged as belonging to one of 
these groups. In order to enhance realism, a 
random variation number for personal data is 
also attached to each individual. The savings in 
storage accrued by this indirect storage assign¬ 
ment are considerable: if we have Nd personal 
data items, conventional storage would require 
O(NdNp) items, whereas indirect data assign¬ 
ment only requires 0(2Np+NdNg) items. This 


implies that even for moderate realism, with 
0(10) data items per person, indirect storage 
provides savings of 0(5). All this data is neces¬ 
sary for more realism in the pedestrian motion 
simulation. 

2.5.2 Geographic data 

The environment where the agents navigate is 
represented in the system as geographic data. 
The geographic data contain information about 
the platform within which the particles move, 
items such as terrain condition (inclination, cor¬ 
ridors, escalators, obstacles, etc.), climate data 
(temperature, humidity, rain, visibility), signs, 
as well as doors, entrances and emergency ex¬ 
its. All this data is stored in a so-called back¬ 
ground grid consisting of triangular elements. 
This background grid is used to define the 
geometry of the problem, and is generated au¬ 
tomatically using the advancing front method 
If40l . All geometrical and climatological data is 
attached to this grid. 

2.5.3 Fixed triangular background grid 

The geography is stored at (fixed) background 
nodes and used to define the geometry. The 
mesh is generated automatically and the data 
is interpolated (linearly) to the pedestrian posi¬ 
tion (see fig. ITT) . 

_Position of Pedestrian 

at Previous Timestep 

! • >r A / 

1 / 'A /’S \ / 

(.—x-A 

• / Position of Pedestrian 

at Current Timestep 

VI_\/ 

Figure 11: Background Grid. All information about the 
terrain surrounding where the particle moves 
will be generated with a fixed triangular back¬ 
ground grid £36}. 

2.6 Pedestrian data 

• Pedestrian "personal" data types 
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- Personal Data (height, width, fit¬ 
ness, nationality, familiarity with 
surroundings, objectives, etc.) 

- Current Physical and Mental State 

- Current Motion Data (location, ve¬ 
locity, etc.) 

• Pedestrian Data Storage 

- Data table for Geno-Types and Ob¬ 
jectives/Paths 

- Pedestrian Arrays for Actual Data 

- Pedestrian "personal" Data types 

2.6.1 Neighbour data 

By far the most time-consuming portion of 
pedestrian simulations is the evaluation of the 
interaction between nearest neighbours. The 
nearest neighbours of every pedestrian must 
be identified and accounted for at every time 
step. Assuming an arbitrary cloud of points, 
a Voronoi tessellation, or its dual, the De- 
launey triangulation, uniquely defines the near¬ 
est neighbours of a point. The Delauney crite¬ 
rion states that the circumcircle of any triangle 
does not contain any other point. Another 
unique way to define a mesh is by minimising 
the maximum angle for any combination of 
triangles adjacent to an edge, (see fig. |12|). 


in parts of the mesh. Every so often (e.g. after 
every time step), the mesh must be modified in 
order to restore it. For 2-D grids like the ones 
contemplated here, this is best achieved by flip¬ 
ping diagonals until the Delauney or min-max 
criterion have been restored. An edge-based 
data structure that is well suited for this pur¬ 
pose stores the two points of the edge, the 
neighbours on either side of the edge, as well 
as the four edges that enclose the edge, (see 
fig. |l~3} . For boundary edges some of these 
items will be missing, making it easy to iden¬ 
tify them. 


Pedestrian 

positions 



Figure 13: Neighbour data. The points represent the 
particle positions and the circles indicate the 
particle-particle interaction ranges. Every 
particle can see and identify the nearest neigh¬ 
bours (at each At). In this system pedestrians 
are moving like a cloud of points, Delauney 
triangulation (updated) BIT . 




Figure 12: Diagonal Swap Using Min-Max Criterion. 

For dynamically moving points, the De¬ 
launey or min-max criterion will be violated 


On the basis of this method the particle can 
identify the nearest neighbours (at each At). 

1.1 Example: walking through a nar¬ 
row passage 

This background grid is used to define the 
geometry of the problem, and is generated au¬ 
tomatically using the advancing front method 
140. 42. 431. All geometrical and climatological 
data are attached to this grid. This implies that 
the amount of data stored and used for pedes¬ 
trian movement depends only on the level of 
details stored in this mesh, and is proportional 
to the number of elements in it. For obvious 
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reasons, the size of this mesh (number of el¬ 
ements and points) should be limited to the 
necessary and available amount of information 
required for the simulation. One can see that 
the mesh will be very fine in areas where the 
density of pedestrians gets higher (see fig. [14) . 
At every instance, a pedestrian will be located 
in one of the elements of the background grid. 
This 'host element' is updated continuously 
with the nearest neighbour tracking procedure. 

Pedestrian motion in passages is one of the 
few cases where reliable empirical data exists. 
In order to assess the validity of the proposed 
pedestrian motion model, a typical 'passage 
flow was selected. The geometry of the prob¬ 


lem is shown in figure 14 Pedestrians enter 
the domain from the left and exit to the right. 
For this case, each pedestrian has the goal of 
reaching first the entrance of the passage, then 
to traverse it to the other end, and finally to 
exit on the right. A typical snapshot during 


one of the simulations is shown in figure 14 





Figure 14: Moving grid and positiori snapshot. This 
picture shoivs pedestrians moving through 
a narrow passage like a cloud of points, De- 
launey triangulation (updated) flit . The fine 
triangulation indicates high density of pedes¬ 
trians. 

The pedestrian parameters chosen were as 
follows: 

• Desired velocity v# = 1 ± 0.02 m/s; 

• Relaxation time r = 0.50 ± 0.1 s; 

• Pedestrian radius R = 0.2 ± 0.02 m. 


The problem was rim repeatedly with an in¬ 
creasing number of pedestrians entering the 
domain per time unit. The resulting pedestrian 
density was measured in the passage, as well 
as the average velocity. The results obtained 


As one can see, the proposed model provides 
surprisingly accurate results. 


are compared with published data in figure 15 



Density (people m 2 /m 2 ) 


Figure 15: Velocity-density fundamental diagram: De¬ 
crease of the pedestrian local velocity v(r, t) 
in the passage area as a function of the local 
density p{f, t). The measured data are repre¬ 
sented by crosses, they are consistent with the 
empirical data of Predtetschenski and Milin- 
ski i(44t represented by the solid line. 


3 Microscopic simulation of 
the Mataf 

The pilgrimage (called Hajj in Arabic) is the 
most significant factor in the life and growth of 
the holy city Mecca. The number of pilgrims 
to Mecca increased in 2009 by 5 percent com¬ 
pared to 2008, official statistics said. In the 
last decade the number of pilgrims has been 
increasing dramatically, every year close to 4 
million pilgrims arrive in the Saudi territory 
through Djedah airport, Djedah port by the 
Red Sea and through the high way connecting 
Saudi Arabia with neighbouring states. The ar¬ 
rival of the huge number of pilgrims each year 
during the last month of the Islamic calendar 
makes this event a grand human gathering in 
the world as well as one of the largest logistical 
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and administrative undertakings [45). The goal 
of pilgrims is to perform the Hajj ritual. Man¬ 
aging the flow of pilgrims to Mecca every year, 
posed various challenges to Saudi authorities 

Si- 

Millions of believers perform the Hajj at the 
same time period, 4 specific days in the year. 
All pilgrims perform the same rituals, typically 
in same or concurrent time periods. The rituals 
are in specific points called points of interest, 
such as the Kaaba in Mecca, the mount Arafat 
etc. From the mountain Arafat, the pilgrims 
continue their course to Muzdalifah to throw 
stones at the pillars symbolizing the devil. Mil¬ 
lions of believers live during the Hajj in tents in 
the Mina area (tent city). After that, millions of 
believers move towards Mecca to do their final 
Tawaf Al-Ifadah and to finish the Hajj ritual. 

In this section we are concerned with Tawaf 
as an important pilgrims activity in Hajj period. 
The accurate prediction of pedestrian motion 
can be used to assess the potential safety haz¬ 
ards and operational performance at events 
where many individuals are gathered. The pre¬ 
diction of the pedestrian flow on the Mataf as 
permanent by overcrowded area and the in¬ 
vestigation of the building facility through the 
simulation can be used to detect the critical 
points with high density in different regions 
of the Mataf. We attempt to determine the av¬ 
erage capacity of the Mataf and we study the 
capacity of the mobile Mataf as possible solu¬ 
tion to reduce the pressure of pedestrian flow 
during the rush hour. As consequence of this 
solution we investigate how the columns of the 
mobile Mataf influence the fluidity of pilgrims 
stream on the ground Mataf. The columns 
serve as a support of the mobile Mataf. An un¬ 
derstanding of how to alter crowd dynamics in 
Hajj activities would have a significant impact 
in a number of other scenarios, e.g. during 
riots or evacuations. Evacuation from danger¬ 
ous zones, restrained places or overcrowded 
buildings, also represent cases where the pre¬ 
diction of pedestrian motion can be used ad¬ 
vantageously. 


3.1 Motivation 

Now that the world population, particularly 
the Muslim population, is increasing, the num¬ 
ber of pilgrims will also increase. With the 
huge development in the transportation tech¬ 
nology in recent years the number of pilgrims 
arriving to Mecca grow systematically. It is a 
fact that the areas of Holy sites, such as Mina, 
Arafat, Muzdalifah and Mecca are fixed. This 
has motivated the Saudi government to find an 
effective and long-term solution. One of these 
solutions is to evaluate and build the mobile 
floor for Mataf expansion. The mobile Mataf 
has a lot of advantages, it can be built in a short 
time on the Mataf area during the Hajj period, 
after that it can be rapidly removed. 

3.2 Observations in the simplified 
model 

One of most important rituals in Hajj is the 
Tawaf, it consists of circling the Kaaba in the 
Holy Mosque seven times counter clockwise. 
Afterwards the Sa'y is performed, the walk¬ 
ing between Safa and Marwa in an enclosed, 
air-conditioned structure. During Hajj seasons 
everything in the Mataf area is dense, and we 
have a compact state. The pilgrims have body 
contact in all directions and they have no influ¬ 
ence on their movement, they are floated by the 
stream. This forms structures and turbulences 
in the flow. This turbulence is very well observ¬ 
able in our video recording f47l . Density and 
velocity can also be seen. The following effects 
could also be observed during the Tawaf: 

• Edge Effects: when the edges of a crowd 
move faster than the center of the crowd. 

• Density Effects: crowd compression in 
local areas can imbalance the crowd flow. 

• Shock Waves: propagation of some ef¬ 
fects spreading throughout a densely 
packed crowd. 

• Speed Effects: higher density causes 
lower walking speeds, with increasing 
contact area between pilgrims. 
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• Group Effects: groups of pilgrims move 
together and try to keep together all of 
the time. 

• Break Out Effects: when pilgrims finish 
the Tawaf, they try to move directly to 
the edge of the Mataf. The movement 
of these pilgrims is normally spiral with 
increasing distance from the center of the 
Kaaba. 

• Structure Effects: With higher density, 
turbulences and structures are created 
within the flow. 

• The average time for Tawaf is ca. 1 hour. 

• While praying in the Haram, there is no 
Tawaf at all. i.e. density decreases when 
praying is going to start and rises when 
praying has finished. 

3.3 The Haram mosque building de¬ 
scription 

The pilgrims stream into the Mataf from all 
gates of the Haram (see fig. [16] (A)), but the 
beginning of Tawaf must start at a specific line 
(see red line in figure [l6|(B)). As we can see in 
figure |16](D), the Haram has 8 gates, stairs and 
an escalator. The gates are probably used by 
all pilgrims. Nearly (approximately 80 percent) 
of pilgrims make Sa'y after Tawaf, the rest (20 
percent) of the pilgrims go out through any 
of the doors. For the simulation all doors are 
treated equally. 

In Tawaf the pilgrims must circle the Kaaba 
seven times in a counter-clockwise direction. 
The black stone designates the beginning and 
end of the Tawaf (the pilgrims must stay 5-7 
s in front of black stone to appreciate saying 
"Bismuallah Allahu Akbar" after every round). 
Before performing Tawaf, pilgrims try to reach 
the black stone. They queue around the Kaaba 
wall with some pilgrims not wanting to queue 
and trying to push through. After the Tawaf 
people try to pray between Maqami Ibrahim 
and the Kaaba (see fig. ??) (The prayer is about 
4 to 5 min). 


In figure [16] (D) a 3D model of the Haram 
Mosque in Mecca is illustrated. This model 
shows the main gates, doors, side entrances, 
and stairs to the Mataf open air area of the 
Haram. The start/end of the Tawaf is indi¬ 
cated by the red line. The blue line indicates 
the Tawaf movement direction on the roof of 
the mosque or in the piazza of the Haram near 
the Kaaba. 

3.4 Simulation of the simplified 
model 

Our first set of simulations consisted of the 
pedestrian flow through the Mataf area with¬ 
out columns (see fig. [17] (B)). The Mataf area 
was 105 m wide and 154 m long and in the cen¬ 
ter of the Mataf the Kaaba with 21.14 m long 
and 11.53 m wide is placed. The circumambu- 
lation area (ground floor) of Mataf is around 
16 170 square meters. Each pedestrian's de¬ 
sired speed was set to a relatively high value 
of 1.2 m/s, the pedestrian radius R was 0.178 
m to 0.2 m and the relaxation time t was set 
in the range of 0.5 s to 1 s. After a pedestrian 
finishes the Tawaf (seven times circling around 
the Kaaba) they go in the direction of Sa'y. The 
measured densities are 5-6 persons/m 2 (mea¬ 
sured in different part of the Mataf area with 
respect to the distance from the Kaaba). These 
results are consistent with those obtained by 
Predtetschenski and Milinski H44l . who also 
found out that high density results in reduced 
velocities. The congested area increases the 
density sometimes up to 7 or 8 persons/m 2 . In 
consequence the pedestrians begin to push to 
increase their personal space and create shock- 
waves propagating through the crowd. 

Figure [17](C) displays a pilgrims movement 
simulation within Mataf area. The entire in¬ 
flux consists of three uni-directional pedestrian 
flows coming from three entrances. The ve¬ 
locity indicator shows that the movement in 
the edges of the Mataf area is faster than in 
the area of the Kaaba. The picture shows a 
snapshot of the simulation at 38 s, which has 
a particularly high density of 6.5 persons/m 2 . 
Also note that the pedestrian density is very 
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Figure 16: (A) Illustration of the Mataf dimension; (B) the direction of the Tawaf movement; (C) the Mataf description: 

(1) the black stone, one of the most overcrowded regions in the Haram, the measured local density p(r, t) 
in this place reached 9 persons/m 2 , (2) the green strips on the ground indicate the begin and end of the 
Tawaf movement, (3) Maqam Ibrahim, hysterical monument near the Kaaba that must be considered in the 
simulation; (D) the gates and entrances to the Mataf. 
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Figure 17: (A) 3D-Model of the mosque building in Mecca. Main gates and entrances are indicated through the red 
arrows, the green arrow indicates the direction of the Tawaf movement, the blue line indicates the beginning 
and the end of the Tawaf and the pink arrow indicates the Sa’y movement; (B) Mataf dimension; (C) PedFlow 
simulation snap-shot results, the numbers 1, 2, 3 indicate the entrances or particles streaming in the Mataf 
area and the number 4 particles streaming out of the Mataf or the exit from the Mataf toward the Sa'y area; 
(D) decrease of the pedestrian velocity in the Mataf area as a function of the local density, the measured data 
are represented by the crosses, they are consistent with the empirical data of Predtetschenski and Milinski 
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high at the places where the Tawaf begins and 
ends. Note the clumping of pedestrians going 
in opposing directions, when the pilgrims fin¬ 
ish the Tawaf. The average density for many 
runs was 5 to 6.5 persons/m 2 , (see fig. 17 (D)). 


3.5 Simulation of the enhanced 
model 

The second set of simulations consisted of the 
pedestrian flow through the Mataf area with 
columns, (see fig. [T8]>. The parameters are 


the same as in section 3.4 After pedestrians 
finish the Tawaf (seven times circling around 
the Kaaba) they are going in the direction of 
Sa'y. The measured densities are again 5-6 
persons/m 2 with up to 7 or 8 persons/m 2 in 
the congested area. Thus, all results are iden¬ 
tical with the simulation without the columns. 


cf. section 3.4 Obviously the columns do not 
influence the densities and velocities of the 
pilgrims in the congested area. The only ob¬ 
servable difference are the vacant rings on the 
circles defined by the positions of the columns. 

Figure [T9] illustrates two simulation snap¬ 
shots with different entire influx. The pedes¬ 
trian streaming the Mataf area from three en¬ 
trances are indicated by the numbers 1, 2, 3 
in figure [l9](B). The velocity indicator shows 
furthermore that the congested area still is in 
the interior zone near the Kaaba and at the 
area where the Tawaf begins and ends. The 
picture shows a snapshot of the simulation at 
38 s, which has a particularly high density of 
6.5 persons/m 2 . 

For the first moment one can not realise the 
difference to the results illustrated in figure [17] 
and figure 18 both pictures exhibit the same 


critical points. Note the clumping of pedes¬ 
trians going in opposing directions, when the 
pilgrims finish the Tawaf. One may conclude 
that the columns have little influence on the 
movement of pilgrims but one can clearly see a 
vacant ring along the columns and this results 
from repulsive interaction between the parti¬ 
cles and the columns, this effect is better seen 
in case of lower crowd density. In case of high 
density crowd the pressure between the indi¬ 


viduals is so huge that this vacant ring along 
the column disappears. This means that under 
the enormous surrounding pressure, people 
are forced to come near the pillar. Here we 
stress that the small size and smooth form of 
the column is very important to prevent people 
injury or dangerous situations. 


4 Mataf capacity estimation 

4.1 Fundamental Diagram 

Pedestrian distribution on the used surface 
must be identified with the fundamental di¬ 
agrams (see fig. |2lJ. 


Velocity v(m/min) 



Different regimes 


Figure 21: Velocity-density diagram; Empirical relation 
between density and velocity according to 
Predtetschenski and Milinski JH- The par¬ 
tition refers to domains with qualitatively 
different decrease of the velocity. The red 
rectangle indicates a dangerous area in this 
curve, where the pedestrian density reaches 
high values with more than 6 persons/m 2 . In 
this range bodies will be pressed or deformed, 
resulting in high injury scales with possible 
fatal consequences. 


4.2 Calculations 

We divided the Mataf area into three parts ac¬ 
cording to the density distribution. The first 
part is the area near the Kaaba wall. The sec- 
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Figure 18: Snapshot 1: (A) 3D Model visualises one favourite concept for developing of the Mobile mataf (see green 
coloured area); (B) CAD drawing displaying the Mataf dimension area with columns; (C) simulation result 
demonstrates how the pillars influence the movement of pilgrims on the Mataf ground. The colour coding 
represents the magnitude of the local velocity, PedFlow snapshot; (D) mobile floor for the expansion of the 
Mataf area in Haram mosque. 
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Figure 19: Snapshot 2: (A) (Mobile Mataf Floor) at Mataf area in Al-Haram Sharif, one can see the columns (pillars) 
bearing the mean structure of the temporary Mataf; (B) A CAD drawing displaying the Mataf area with 
columns; (C) and (D) illustrating a movement of pilgrims in the Mataf area (PedFlow simulation results). 
The simulation shoivs a vacant ring since the pilgrims avoid the circle on which the columns are placed. In 
the region of the beginning and end of the Tawaf ingoing people push back the outgoing people. Thus, it 
creates a trap becoming more and more crowded. 
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Figure 20: Pilgrims streaming the Mataf area through three different gates denoted by 1, 2, 3 and leaving the Mataf 
floor through gate 4. Step by step simulation of a pedestrian flow circling seven time the holy Kaaba ((A), 
(B), (C) and (D) illustrate time consecutive PedFlozv animation snapshots). The impact of the columns on 
the pilgrims movement performing the Tawaf are demonstrate. The red colour indicates pedestrian walking 
with their desired velocities and the turquoise colour indicates pedestrian walking with lower velocities. 
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ond part is the middle of the Mataf, the area 
between 10 m and 20 m from the Kaaba wall. 
The third part is the edge of the Mataf, the area 
between 20 m and 30 m. 

• The first part of the Mataf area: 

- Usable Area = 820 m 2 . 

- The mean pedestrian density = 6 
pilgrims/m 2 . 

- The total number of pilgrims = 4920 
pilgrims. 

- The mean pedestrian path = 91 m. 

- The mean Tawaf path = 91 x 7 = 637 

m. 

- The Tawaf lasts 39 min with aver¬ 
age walking velocity of 0.266 m/s = 
15,96 m/min. 

- The capacity is 7596 pilgrims 
Tawaf/hour 

• The second part of the Mataf area 

- Usable Area = 1380 m 2 . 

- The mean pedestrian density = 3,68 
pilgrims/m 2 . 

- The total number of pilgrims = 5078 
pilgrims. 

- The mean pedestrian path = 151 m. 

- The mean Tawaf path = 151 x 7 = 
1057 m. 

- The Tawaf takes 52,85 min with av¬ 
erage walking velocity of 0,33 m/s 
= 20 m/min. 

- The capacity is 5859 pilgrims 
Tawaf/hour. 

• The third part of the Mataf area 

- Usable Area = 1900 m 2 . 

- The mean pedestrian density = 2 
pilgrims/m 2 . 

- The total number of pedestrians = 
3800 pilgrims. 

- The mean pedestrian path = 212 m. 


- The mean Tawaf path = 212 x 7 = 
1484 m. 

- The Tawaf lasts 49,46 min with av¬ 
erage walking velocity of 0.5 m/s = 
30 m/min. 

- The capacity is 4653 pilgrims 
Tawaf/hour. 

The total capacity of the Mataf is 18923 pil¬ 
grims Tawaf/hour. 

According to several studies about the 
Haram, the actual capacity of the entire Mataf 
(ground, intermediate floors and the roof) are 
around 50,000 pilgrims/hour, but the capacity 
of the new Massaa is much higher than the 
capacity of the existing Mataf. During peak 
times worshippers are in danger of injuries 
due to extreme high density in the Mataf. To 
increase the capacity and improve the condi¬ 
tions for worshippers a new additional floor 
for the Mataf called temporary Mataf has been 
developed. 

5 Mobile Mataf geometry 

The new floor is based on a modules concept 
with fast assembly and disassembly without in¬ 
terruption of ongoing Tawaf. This light weight 
structure is designed to be integrated into the 
existing Haram Mosque. The outer diameter is 
about 94 m, the inside diameter is around 70 m 
and the width is around 12 m (see fig. [23] (A)). 
An additional area of 3000 m 2 is added to the 
existing Mataf area. Access is achieved over 
several ramps. The clearance height is around 
2.7 m leaving enough space for worshipper to 
pass underneath. 

5.1 Concepts and evaluations 

The idea to develop a new moving platform 
above the existing Mataf has the goal to take 
pedestrian pressure from the old Mataf. The de¬ 
sign and evolution of this mobile area revealed 
a new concept concerning the safety and flu¬ 
idity of pilgrims flow. In the next paragraphs 
we discuss two approaches with two different 
geometries and we present our analysis data. 


19 



Running title • November 2012 • Vol. XXI, No. 1 


• Mobile Mataf with extension: we shall 
discuss 


- Individual Radius: R = 0.178 ± 0.02 
m (Milinski). 


- Geometry 

- Simulation results 

- Capacity of the Mataf according to 
Predtetschenski and Milinski densi¬ 
ties 

5.2 Simulation 


- Relaxation Time: r = 0.5 ± 0.1 s. 
Geometry Parameters. 

- Usable Area: A = 3540 m 2 . 

- The mean Radius of the Mataf: R 
= 42 m. 


This set of simulations considers a unidirec¬ 
tional pedestrian flow circling the Holy Kaaba 
seven times. Different strengths of the influx 
of pilgrims were used, see table [2]The walking 
pedestrian area named mobile floor for the ex¬ 
pansion of the Mataf contains two exits and 
one emergency exit. The design of the Mo¬ 
bile Mataf allows one-directional pedestrian 
flow (see fig. 23 (A)). The parameters used 
in the simulations were the desired velocity 
of the individual V&, the individual radius R 
and the relaxation time T. The size or radius 
of particles are given in Table |T] and accord¬ 
ing to Predtetschenski and Milinski l44l . For 
simplicity, all parameters are identical for each 
person. 


Type of Person 

Projected area 
(m 2 ) 

Children 

0.04 -0.06 

Adolescent 

0.06 -0.09 

Grown-up in summer 
clothes 

0.100 

Grown-up in interseason 

0.113 

Grown-up in winter clothes 

0.125 

Grown-up in interseason 
clothes with briefcase 

0.180 

Grown-up in interseason 
clothes with light luggage 

0.240 

Grown-up in interseason 
clothes with heavy luggage 

0.390 


Table 1: Set of different individual sizes Hit . 

5.2.1 Parameters 

• Individual Parameters. 

- Desired velocity: Vj = 1.2 ± 0.2 m/s. 


- The mean pedestrian path: L Pat f, = 
2 tcRm = 263.76 m. 

- The mean Tawaf path: Lj awa f = 
263.7 x 7 = 1846.32 m. 

- The Tawaf time: tj awa f = 43.35 min 
with average walking velocity v = 
0.7 m/s = 42 m/min. 


5.2.2 Results 


Influx (Per¬ 
sons/s) 

Max Num¬ 
ber of Ped 

Density (1/m^) 

Capacity (Tawaf 

h) Remarques 

2,0 

5032 

1.42 

7021 

fluid 

2,1 

5338 

1.50 

7448 

fluid 

2,2 

5796 

1.63 

8087 

fluid 

2,3 

6147 

1.73 

8577 

fluid 

2,4 

6691 

1.89 

9336 

fluid 

2,5 

6955 

1.96 

9704 

fluid 

2,6 

7443 

2.10 

10385 

fluid 

2,7 

8106 

2.28 

11310 

fluid 

2,8 

8519 

2.40 

11886 

fluid 

3,9 

9236 

2.60 

12887 

Slow Move¬ 
ment 

3,0 

10358 

2.92 

14453 

Clogging ef¬ 
fect 


Table 2: The Mataf capacity at different simulation in¬ 
fluxes. 


In table [2] the maximum number of pedestrians 
on the Mataf without clogging is 8404 with a 
density of 2.37 persons/m 2 . The optimal ca¬ 
pacity without clogging is 11726 pilgrims per 
hour completing the Tawaf, calculated with an 
average Tawaf duration of 43 min. 

The maximum number of pedestrians on 
the Mataf was 11629 with a density of 2.79 
persons/m 2 . The optimal capacity without 
clogging was 16226 Tawaf/hour, calculated 
with an average Tawaf duration of 43 minutes. 
If one Tawaf is 50 minutes then the capacity of 
the Mataf will be 13954 Tawaf/hour. 
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Figure 22: (A) CAD drawing displaying one of many geometry concepts developed for evaluation of the mobile Matafi 
(B) simulation results (new geometry concept); (C) velocity density diagram (PedFlow simulation); (D) the 
actual mobile floor facility top view. 
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6 Mataf capacities estimation 


According to the Predtetschenski and Milinski 


fundamental diagram figure 21 the capacity 
of the mobile Mataf was determined. The re¬ 
sults of this investigation are illustrated in Ta¬ 
ble [3] The capacities of the mobile floor were 
computed with respect to different pedestrian 
densities. Every density in the fundamental 
diagram is related to a certain velocity. With 
help of these given velocities we tried to calcu¬ 
late the mean average time for one complete 
Tawaf. Afterwards we calculated the capacity 
of the Mataf. We can see that the optimal ca¬ 
pacity is about 12000 Tawaf/hour, (see fig. 23 
(C)). This result consistent with the simulation 
result which is about 13954 pilgrims per hour 
completing the Tawaf. 


Density 

(1/m 2 ) 

Velocities 

(m/min) 

Path 

length 

(m) 

Tawaf 

time 

(min) 

Usable 

Area 

(m 2 ) 

Max 

Pedes¬ 

trian 

num¬ 

ber 

Capacity 
(Tawaf/h 

0 

60 

1846 


3864 

0 

0 

0.05 

47.19 

1846 

39.11 

3864 

1932 

2963 

0.1 

38 

1846 

48.57 

3864 

3864 

4773 

0.15 

32.98 

1846 

55.97 

3864 

5796 

6213 

0.2 

30 

1846 

61.53 

3864 

7728 

7535 

0.3 

22 

1846 

83.90 

3864 

11592 

8289 

0.4 

18.19 

1846 

101.48 

3864 

15456 

9138 

0.72 

13.5 

1846 

136.74 

3864 

27820 

12207 

0.92 

9.03 

1846 

200.65 

3864 

35548 

10629 


Mataf area. The new Mataf area was designed 
to satisfy the international standard criterion 
of pedestrians safety, especially during over¬ 
crowding. The flow of pilgrims is demanding, 
since any accumulation in this area can be dan¬ 
gerous and has fatal consequences. The new 
proposed design offers a significant improve¬ 
ment to the safety of pilgrims. The limited 
spatial possibilities in and around the Haram 
area for infrastructural development requires 
an in-depth analysis and search for solutions 
to solve the problem of pedestrian flow. From 
this the concept arises to build a temporary 
Mataf to take a small capacity from the ground 
Mataf as best arrangements during the massive 
extension project in Masjid Al-Haram. 



Table 3: Mataf capacities according to Predtetschenski 
and Milinski densities. 


Figure 24: The new walk way around the Kaaba 
is definitely taking shape in this picture 
(26.07.2013). 


7 Recommendations for the 
future Hajj planning and 
development 

With the increasing number of pilgrims every 
year and with the dramatic evolution of the 
transportation system in the world, the de¬ 
sign and architecture of the Haram mosque 
in Mecca must be changed to increase the 
capacity, throughput and performance. The 
proposed and evaluated design has a specific 
aim, the safety of pilgrims. In order to en¬ 
sure that the proposed design satisfies the re¬ 
quired criteria of pilgrims safety, especially dur¬ 
ing overcrowding, many simulations and tests 
were realized to access the critical points in the 


The Mataf area is a very restrained plain. 
The construction of the temporary Mataf 
within Masjid A1 Haram should attempt to 
resolve many problems to avoid a potential dis¬ 
aster in the future, such as the overcrowding 
disaster at the Jamarat bridge in 2005. The con¬ 
ceptual design as we can see in figure |24| has 
two problems; the first problem concerns the 
fluidity of the movement on the platform espe¬ 
cially at the exits, the second problem concerns 
the ground Mataf in that the fluidity of pedes¬ 
trian flow will be disturbed through the pillars 
carrying the platform. The Haram building 
must ensure the international norm of safety 
to allow the increased number of pilgrims. 

To prevent stampede and congestion within 
the Mataf area during the rush hour of the 


22 



























Running title • November 2012 • Vol. XXI, No. 1 



Figure 23: (A) CAD drawing displaying the geometry concept of the new walkway around the Kaaba; (B) Predtetschen- 
ski and Milinski fundamental diagram; (C) mobile floor capacity according to Predtetschenski and Milinski 
densities; (D) the present mobile floor facility top view. 
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Hajj events a sufficient study and analysis of 
the expected inflows and outflows (and, hence, 
number of pilgrims) must be ascertained and 
documented, considering the alternation possi¬ 
bility of large pedestrian flows. A narrow pas¬ 
sage and bottleneck analysis within the Haram 
mosque area is crucial for any decision mak¬ 
ing. Through exact simulation and analysis of 
pedestrian fluidity on the Haram area the con¬ 
gestion points as the intersection or crossing 
points must be established and should be re¬ 
moved. Either the passages with counter flow 
must be detected as early as possible and care¬ 
fully treated (examined). This analysis can be 
used to assess the potential safety hazards and 
can be helpful for the organizers. As critical 
points the exits of the Haram must be checked. 
The passage between the Mataf area and the 
Massaa area and between the pillars of the 
temporary Mataf can cause problems for the 
fluidity of the pilgrim streams. 

For the safety of all pilgrims on the Haram 
area and to ensure a quick and effective re¬ 
sponse to any occurring problems a sophisti¬ 
cated and continuously monitoring system (e.g. 
video surveillance) must be established. 

An engineering system must be developed 
within the Haram to allow the security and 
emergency forces to move as quickly as possi¬ 
ble and unobstructed to any place if something 
happens, to remove or at least mitigate prob¬ 
lems. 

The main problem in Mecca itself are the 
traffic jams. The high density of cars and buses 
combined with walking or road crossing pedes¬ 
trians make the situation more complex for res¬ 
cue or evacuation plans. The possibility to be 
stampeded at Hajj pilgrimage is very high at 
all places and especially inside the Haram in 
Mecca, in Muzdalifa, at the Jamarat and on the 
streets between the so called points of interest. 

The changes necessary for a successful Hajj 
and its activities are all contained in the wider 
long-term plan for Mecca. The next expan¬ 
sion for Haram and Massaa requires a new 
infrastructure and regeneration of the east end 
- which will reclaim old neighbourhoods to 
Haram for housing and parking area - which 


is a huge project. 

8 Conclusions and possible 
enhancements 

A newly developed simulation software for 
analysing pilgrim movement was used for the 
evaluation of the Tawaf capacity of the Mataf 
area and the new mobile floor. Different sce¬ 
narios were analysed to find out the optimal 
capacity. With a flux of only 12 pilgrim/min 
entering the floor, much space was left un¬ 
used. Pilgrims moved with their desired speed 
and the average time for completing a Tawaf 
took around 25 minutes. The second scenario 
showed 108 pilgrim/min entering the floor, 
causing the space to become fully occupied 
by pilgrims and the velocity to decrease, re¬ 
sulting in dangerous jams especially in the exit 
area. The third scenario showed the optimum 
steady state when about 125 pilgrim enter per 
minute and also exit the floor, with no jams. 
This results in an overall capacity of about 9000 
finished Tawaf/hour and an average time of 
approximately 25 min to complete the Tawaf. 

This work is an attempt to model high den¬ 
sity crowd dynamic flows in autonomous and 
controlled scenarios. A microscopic approach 
was used to model this problem based on ob¬ 
served collective behaviour in emergency con¬ 
ditions, where the detailed design of interac¬ 
tions is overlapped by group behaviour. We 
presented a simple means of achieving more 
complex behaviour in both indoor and outdoor 
environments. The simulation tools produce 
results that compare favourably with the real 
data. 

As a practical example, the Haram Mosque 
in Mecca and the Jamarat Bridge in Saudi Ara¬ 
bia were used for high density crowd simula¬ 
tion: the huge number of pilgrims cramming 
the bridge during the pilgrimage to Mecca 
caused serious pedestrian disasters in the 90's. 
Moreover, the analytical and numerical study 
of the qualitative behaviour of human individ¬ 
uals in a crowd with high densities can also 
improve traditional socio-biological investiga- 
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tion methods. 

For obtaining empirical data different meth¬ 
ods were used, automatic and manual methods. 
We analysed video recordings of the crowd 
movement in the Tawaf within the mosque in 
Mecca during the Hajj of the 27th of November, 
2009. We evaluated unique video recordings of 
a 105 x 154 m large Mataf area taken from the 
roof of the Mosque |47| . 

For the validation and calibration of the 
simulation tools, different methods were used 
S48J. At medium to high pedestrian densities, 
the techniques used in PedFlow can produce re¬ 
alistic crowd motion, with pedestrians moving 
at different speeds and under different circum¬ 
stances, following believable trails and taking 
sensible avoidance action. 

As possible future work, we would like to 
have a better understanding of the parameters 
influencing the crowd behaviour such as the 
ethnic and psychological parameters. These 
parameters have a significant impact on the 
movement of crowds especially in panic sit¬ 
uations. The psychological-physiological pa¬ 
rameters in the case of riots or emergencies 
must be investigated further and documented 
so that they can be implemented into other 
models. Our proposed approach to simulate 
crowd behaviours in the future is based on this 
dichotomy, and we propose separate models 
for an agent's personality and another one to 
account for situational factors. This question 
arises and opens the door to many topics that 
must be further investigated. There are certain 
things we can do to help a high density crowd 
participant like pilgrims, sport or music event 
spectators to survive a stressful situation, by 
planning, preparation, training and better man¬ 
agement. But this is not sufficient to make sure 
that hundred percent of the crowd are safe, if 
we consider the number of accidents happened 
in the last decade. 
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